The time dependence of muon spin relaxation spectra and spin correlation functions 
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The existing theory of the microscopic interpretation of the dynamical contribution to zero-field 
muon depolarization spectra in a longitudinal geometry is developed. The predicted relaxation of 
the muon depolarization is calculated from two forms of the spin correlation function. First, when 
the spin correlation function has an exponential form with a single wave vector dependent relaxation 
rate is considered, it is shown that this form of the spin correlation function reproduces the slow and 
fast fluctuation limits of stochastic spin theory regardless of the choice of microscopic spin model. 
Second, if the spin correlation function is a homogeneous scaling function (such as a power-law decay 
with time), as suggested by the mode-coupling theory of spin dynamics, this results in a stretched 
exponential relaxation of the muon spectra. For simple spin diffusion, the muon spectra are shown 
to be relax with a root-exponential form. 
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I. INTRODUCTION 

The study of spin dynamics is fundamental to the clas- 
sification of magnetic materials. One tool that has been 
used extensively is neutron scattering which has pro- 
vided experimentalists with detailed knowledge of the 
both static and dynamic properties of magnetic systems. 
However neutron scattering is inappropriate for samples 
with a large absorption cross-section, or which are only 
available in small quantities, or when the magnetic mo- 
ment is small. 

Muon spin relaxation spectroscopy (/iSR) is an com- 
plementary probe of the spin dynamics in magnetic ma- 
terials, and can be used for materials when neutron spec- 
troscopy is not suitable^. The muon is implanted into a 
magnetic lattice and is sensitive to the spin fluctuations 
through the fluctuations of the local magnetic field. The 
motion of the muon spin is sensed when it decays into 
a positron and other products. The positron is emit- 
ted preferentially in the direction of the muon spin at 
the time of its decay. This can be studied in longitu- 
dinal muon relaxation spectroscopy where the emitted 
positrons are counted in two banks of detectors, one in 
front and one behind the sample position. The quan- 
tity of experimental interest is the time-dependent muon 
asymmetry which measures the time dependence of the 
motion of the muon depolarization. This is given by 



Pz{t) 



N F (t)-aN B (t) 
N F (t)+aN B (t)' 



(1) 



where Np(t) is the number of positron counts in front 
bank of detectors, Nsit) is the number of positron counts 
in the rear bank of detectors and a is an experimentally 
determined constant. 

The longitudinal asymmetry function, Pz(t), can have 
static and dynamic contributions. In the case of com- 
pletely static local fields, the muon senses an average of 
the local field given by 



Pz(t)= J p(B^{cos 2 6 + sm 2 (6)cos{ lf ,B t ,t)}dB f „ 



(2) 



where 9 is the angle between the local field (given by the 
vector B^) and the initial muon spin, p(B M ) is the mag- 
netic field distribution function and 7 M is the muon's gy- 
romagnetic ratio (2irx 13.55 kHz /G). The choice of field 
distributions involves some motivation from the knowl- 
edge of the magnetic system. For instance for a dilute 
spin glass the field distribution is likely to be a Gaus- 
sian distribution, which gives rise to a Kubo-Toyabe 
relaxation function^. One other relaxation function of 
note, as it occurs often, is the stretched exponential 3 
(Pz(t) cx exp {— (At) K }). The form where k = \ is often 
seen in the low temperature properties of various types of 
magnetic material a 4 i 5l6 i and is attributed to the Uemura 
distribution of local fields^. 

This paper is concerned with the amount of informa- 
tion that can be obtained about dynamic spin fluctua- 
tions from muon spectra. Previous theoretical studies of 
the dynamic contribution to the shape of the muon asym- 
metry have predicted a simple exponential form, as they 
have been studied in the fast fluctuation limit (when the 
spins fluctuate over times much shorter than the muon 
experimental time window) 8 . These studies established 
a direct relationship of the spin correlation function and 
the muon asymmetry However, it has been known for 
some time that the spin correlation function may exhibit 
long-time power-law tails&iS that may be significant in 
the muon time window. 

It is possible to identify two regimes for spin cor- 
relation function, the high temperature paramagnetic 
regime, where the spin correlation function exponentially 
decaysii and the critical regime where it is a homoge- 
neous function^ of the scaling variable q~ e t (where q is 
the wave vector, t is time and is an exponent) 10 . In this 
paper the dynamic contribution to the time-dependence 
of the muon asymmetry is derived, without using the fast 
fluctuation limit, for these two common forms of the spin 
correlation function. 

In the next section the theory describing how spin fluc- 
tuations affect muon spectra are briefly reviewed. A 
method of Laplace transforming the muon polarization 
function is used and this is applied to the case when the 
spin fluctuations are decaying exponentially with time. 
Mode-coupling theory suggests that the spin correlation 
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function should be a homogeneous function of the prod- 
uct q t (where q is the wave vector and t is time) so 
the effect on the depolarization spectra of this form is 
studied. It is shown that this can result in a stretched 
exponential relaxation and, in the case of simple spin dif- 
fusion, root exponential behavior. The exponent (k) is 
dependent on the model of spin dynamics chosen. 



II. A BRIEF REVIEW OF EXISTING THEORY 

The muon depolarization shape has been calculated 
using perturbation theory. A quantum calculation 13 and 
a semi-classical calculation 14 lead to 



P z (t) =exp 



(3) 



where the coordinate axes are such that the z-axis is 
aligned with the initial muon spin direction. The relax- 
ation function \t z is given by 



f dr{t-T)[$ a 

Jo 



(r) + %y(r)}, (4) 



where $ aa (t) = \{5B a (t)5B a + SB a SB a (t)) is the field 
fluctuation correlation function for component a of the 
local magnetic field B a . This is a result of perturbation 
theory so it is valid when the interaction term of the 
muon and local field is much smaller than the other terms 
in the Hamiltonian. 

The dynamical behavior of the local field has been 
analysed before in the fast fluctuation limit when the field 
correlation functions are significant over a much shorter 
timescale than the muon experimental time window 8 . 
This means that one can say that equation 0] is propor- 
tional to time and the relaxation rate (A z ), in this fast 
fluctuation limit, is given by£ 



A 2 



d 3 q 



J 01 v 



(5) 



where V = (/io7^3A i s/47r) 2 , A is a tensor describing how 
the spins couple to the local magnetic field, A(q, w M = 0) 
is the q-dependent spin correlation function approxi- 
mately at zero-energy because of the low energy nature 
of the muon probe, g is the usual Lande g- factor, /is is 
the Bohr magneton, 7^ is the muon gyromagnetic ratio 
and V is the volume of the sample. Also the wave vector 
integration in equation [S] is over the first Brillouin zone. 

The time-dependence of the muon depolarization func- 
tion, P z {t), is not necessarily just a simple exponential 
and this paper is concerned with the predictions of the 
shape of the muon depolarization equation that can be 
made from equation Laplace transforming equation 0] 
gives 

*.(*)= 7^ {***(*) + * w 0>)}- (6) 
where s is the Laplace conjugate variable to t. The 
Laplace-transformed field correlation functions are then 



given by 



V 



{<Ms) + <M*)} = {^} 21 



(7) 



P,1=z,V,z 



where A?"* = G xf3 (q)G x ~< (-q) + G y/3 {q)G v ^ (-q) (where 
(3 = x,y,z and 7 = x,y,z) and D is the dimensional- 
ity. This function G a @ is the Fourier transform of the 
dipole-dipole and hyperfine interactions between muon 
and lattice spins. Near q ~ this is given by 8 



G a P(q) = -47T \ Pf(q) - C Q/3 (0) - ^ 



(8) 



where H is the Hyperfine field at the muon site, r M is 
the number of nearest neighbours to the muon, C a "(0) 
is a symmetric tensos^ and P"^(q) — q a qP /q 2 . In this 
long wavelength limit the coupling tensor A& 1 is solely 
dependent on the angular variables of the integral. Using 
spherical coordinates in the integral gives 



fQBZ 



Mo \ 2 1 
4^W V 



E 



(9) 



where the angular part of the integration is contained 
within the co-efficient B@~ f , which is given by 



g/37 = / A^dttq. 



(10) 



III. 



EXPONENTIALLY DECAYING SPIN 
CORRELATIONS 



At high temperatures the spin correlation function for 
a simple ferromagnet has been shown to be an exponen- 
tial function of tim o 9 ! 11 . Therefore, the spin correlation 
function is given by 



A^(q,t) = 



Nk B T 
9 Mo/4 



0i 



(11) 



where N is the number of spins, T is the temperature 
and r' 37 (q) is the (/-dependent relaxation rate. Laplace 
transforming and substituting this into equations El and 
El gives 



2 4+D 7: 2+D v E 

0,1=x,y,z 



ffPl 



T^(q)} 



dq. 



(12) 



where v is the volume per spin. Inverse Laplace trans- 
forming this gives 



3 



0,1=x,y,z 



qBZ ^-y^) 
[r^(?)] 2 

i 



1 (<?)* 



•r^(g)t}dg. 



(13) 



where = ( 7 2 / i fc B T)/(2 4 — tt— 

muon relaxation function is a q-integrated Abragam 
function. 1 ?. When the relaxation rate is independent 
of wave vector, this result also reduces to the result of 
KereniS in zero-field. The present extension however de- 
rives this result from a microscopic viewpoint as sug- 
gested by Yaouanc and Dalmas de Roetier— . 

Stochastic spin dynamical theoryAS predicts that in the 
fast and slow fluctuation limits the muon depolarization 
spectra are a simple exponential and Gaussian respec- 
tively. It is possible to consider the fast and slow fluctu- 
ation limits of equation 1131 In the fast fluctuation limit, 
which is when T^ 1 {q)t 3> 1, the expression in the curly 
brackets of equation 1131 tends to T(q)t so the relaxation 
function reduces to fast fluctuation limit expression^, and 
is given by 



The derived 



0,1 



-i X^(g) 



dq. (14) 



In this limit the depolarization spectra are a simple ex- 
ponential. In the slow fluctuation limit (T^ 7 (q)t -C 1) 
the expression in the curly brackets of ea uat ion 1 1 31 tends 
to h(T(q)t) 2 so the relaxation function becomes 



*.(*) 



X^(q)q n - 1 dq} . (15) 



This analysis suggests that the spectra can be analysed 
continuously using equation 1131 as the fluctuations slow 
down and pass through the muon experimental time win- 
dow. The transition between these two limits will be 
dependent on the spin model chosen. 

In order to gain some insight the contribution to the 
muon relaxation spectra from one component can be cal- 
culated. The depolarization function can be written as 

P z (t)=P z (Q)]Je X p{-^(t)}. (16) 

0i 

Here \f 7 (t) is given by equation [T3l for one value of ,$7. 
Then we use the Ornstein-Zernike form of the suscepti- 
bility x Pl (q) = X /37 (0)C 2 /( 1 + (<?C) 2 ), where £ is the mag- 
netic correlation length. This form for the g-dependent 
susceptibility has been used in the derivation of equation 
[SI and in previous studies of spin dynamics with /iSR 8 . 

Also for three dimensional paramagnets at high tem- 
peratures, Blume and HubbarcUi found that the spin cor- 
relation function, at long times, has a fluctuation rate 
given by T /3 ' y (q) = Tq 2 (where T is a fundamental spin 
fluctuation rate). The prefactors (^i B^ 7 x^7(o)) are set 




rq= 7 (ns-') 



FIG. 1: The contribution to the muon relaxation function 
from one spin component for a three-dimensional paramag- 
net obeying T' 37 (g) = Tq 2 . The transition from a Gaussian 
to simple exponential relaxation function can be seen as the 
scaled spin fluctuation rate (rq% z ) increases. This is plotted 
for a single value of the magnetic correlation length q£ — 0.1. 



to 1 as we are only interested in the shape change as the 
system passes through a spin-freezing transition. 

Using this and numerically integrating over the Bril- 
louin zone in eouation ll3l we can find the normalised re- 
laxation function exp 

{-*07(t)} which is plotted in fig- 
urc^for one value of the correlation length (qBz£, = 0.1). 
It can be seen that for high values of the spin fluctu- 
ation rate the relaxation is flat and simple exponential 
and that there is a transition to Gaussian behavior be- 
low Tq 2 BZ ~ 10//S -1 . 



IV. POWER LAW DECAY OF SPIN 
CORRELATIONS 

The mode-coupling theory of spin dynamics provides 
an accurate description of the long wavelength proper- 
ties of spin fluctuation9iS*A2i2L2i. The spin correlation 
function is predicted to be a homogeneous function of 
the variable q t for long wavelengths and times, where 9 
depends on the spin model chosen. This scaling has also 
been predicted for a simple ferromagnet near criticality 9 
for 9 = 5/2. Using the spin correlation function given by 



A^(q,t) 



k B TN 



X 01 {q)f^{q% (17) 
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Then the Laplace transform of this is 

A ^ {qfS) = J^L x ^ {q)q -e^ {q -e s)) 
9 MoMs 



(18) 



where q~ e J-^ 1 (q~ e s) is the Laplace transform of 
f^ 7 (q e t). This is a general form of the spin correla- 
tion function and does not include the exponential form 
of the preceeding section as that is not a homogeneous 
function- 2 . It is, however, a general form of the power- 
law time decay of the spin correlation function that is 
predicted for simple spin diffusioniS. The Laplace trans- 
formed shape function is given by 



*,(*) = *i 0) 



fry X^(0) 
Q s 3-(D/8) 



i + e 2 (f)' 



-dx, 



(19) 



where x = q~ e s. The integral, assuming it is finite and 
exists, for small s becomes independent of s. In this limit 
the shape function tends to 



r/37 

$ ( s ) ~ $(o) V B Pl -2— 



where the integral has been replaced by J^- 7 . Equation 
E01 has the simple inverse transform of 



(20) 



07 



9M) ~ $(0)t2-(O/fl) y £07 _ fZ 

,™ er(3-(D/S)) 



, (21) 



where r(x) is the Gamma function and not a relaxation 
rate. 

Spin diffusioniS gives a value of 9 = 2, which in three 
dimensions gives root exponential relaxation 



P z (t) ~P z (0)exp{-A 1/2 Vt}, 



(22) 



where Ajya is given by the prefactor of t 2- ^/ 6 ) in equa- 
tion for D = 3 and = 2. This root exponential 
relaxation has not been predicted before from a dynam- 
ical approach but has been observed many times^<£<£. 

For a three-dimensional ferromagnet near its critical 
temperature, 9 = 5/2 and D = 3, this predicts that the 
muon will relax with the time-dependence 



P z (7i)~P z (0)exp{-A4/ 5 t 4/5 } 



(23) 



It would be an interesting test if this form of the muon 
depolarization spectra could be observed in a model fer- 
romagnet. 

Therefore this suggests that when stretched expo- 
nential relaxation of the muon polarization is observed, 



the spin correlation function is likely to be of the form 
of equation 1171 and possibly a power-law time decay. 
The stretched exponent is also directly related to the 
microscopic form of the spin correlation function. 

V. CONCLUSIONS 

A method that develops the theory of spin fluctuations 
and their effect on the dynamic contribution to muon 
asymmetry has been presented. Two forms of the spin 
correlation function have been considered, corresponding 
to high temperatures in the paramagnetic regime and 
when the magnetic system is in a scaling regime (close to 
criticality). 

In the first case the spin correlation function was con- 
sidered to be a simple exponential with a g-dependent 
relaxation rate. The muon depolarization spectrum was 
found to be an g-integrated Abragam function (equation 
II 3|) that in the fast and slow fluctuation limits recovered 
the forms predicted by stochastic spin theory 18 . The 
asymmetry spectra were seen to be a simple exponential 
when the spin fluctuation rate is large, say at high tem- 
peratures in the paramagnetic regime, and exhibited a 
fast drop when the shape changed to a Gaussian when 
the spin fluctuation rate slowed down. The transition 
between these two limits depends on the specifics of the 
spin model chosen. 

The second case considered was one where the spin cor- 
relation function is a homogeneous function of a scaling 
variable q e t. This form was suggested by the spin diffu- 
sion and mode-coupling theories and is consistent with a 
power-law time-decay of the spin correlation function. In 
this case the asymmetry is found to be a stretched expo- 
nential with a stretched exponent of k = 2 — (D/9) where 
D is the dimensionality and 8 is a coefficient that depends 
on the microscopic model of the spins. Root exponential 
relaxation occurs in three dimensions for simple spin dif- 
fusion and k = 4/5 for a three-dimensional ferromagnet 
in its scaling regime. 

Muons provide a useful alternative probe of spin fluctu- 
ations. The data can often be obtained quicker than with 
neutrons and for samples where neutron spectroscopy is 
not suitable. It is hoped that it has been demonstrated 
that the shape of the muon depolarization spectrum con- 
tains useful information about the type of spin correla- 
tions present and that the analysis presented will help to 
stimulate discussion in this area. 
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